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ABSTRACT 

We present detailed Fourier modeling of the radio remnant of Supernova 
1987A, using high- resolution 9 GHz and 18 GHz data taken with the Australia 
Telescope Compact Array over the period 1992 to 2008. We develop a parame- 
terized three-dimensional torus model for the expanding radio shell, in which the 
emission is confined to an inclined equatorial belt; our model also incorporates 
both a correction for light travel-time effects and an overall east-west gradient in 
the radio emissivity. By deriving an analytic expression for the two-dimensional 
Fourier transform of the projected three-dimensional brightness distribution, we 
can fit our spatial model directly to the interferometric visibility data. This 
provides robust estimates to the radio morphology at each epoch. The best-fit 
results suggest a constant remnant expansion at 4000 ± 400 km s -1 over the 16- 
year period covered by the observations. The model fits also indicate substantial 
mid-latitude emission, extending to ±40° on either side of the equatorial plane. 
This likely corresponds to the extra-planar structure seen in Ha and Lya emission 
from the supernova reverse shock, and broadly supports hydrodynamic models in 
which the complex circumstellar environment was produced by a progression of 
interacting winds from the progenitor. Our model quantifies the clear asymmetry 
seen in the radio images: we find that the eastern half of the radio remnant is 
consistently ~ 40% brighter than the western half at all epochs, which may result 
from an asymmetry in the ejecta distribution between these two hemispheres. 
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Introduction 



Supernova (SN) 1987A in the Large Magellanic Cloud (at a distance of 51 kpc) was 
the nearest supernova observed in almost 400 years, and has enabled detailed studies of 
supernova physics. The detection of neutrino bursts indicated that the event was a core 
collapse, however SN 1987A was not a typical example of a type II supernova in many aspects- 



Optic al observations revealed a triple-ring nebula surrounding the explosion (IBurrows et al. 



19951 ). This indicates the inner surface of the densest circumstellar material (CSM) that 
was ionized by the supernova UV flash. Observations of the l ight echos provided a more 
comprehensive picture of the structure (ISugerman et al.l l2005al Jbl) . The CSM was found to 
have a peanut-like bipolar geometry, inclined at 43° to the line of sight. 

The complex circumstellar environment can be understood as analogous to t he in- 
terac ting stellar wind model, as originally proposed to explain planetary nebulae (IKwok 
19821 ). The progenitor star of SN 1987A, Sk — 69°202, is believed to have evolved from 
a red supergiant (RSG) to a bl ue supergiant (BSG) about 20000 years before core col- 
lapse (ICrotts fe Heathcotd Il99ll ). After the transition, the fast, low density BSG wind 
overran the slow, dense RSG wind and swept it into a circumstellar shell. Some mod- 
els suggest that the RSG wind could have been highly aspherical, due to a binary m erger 
( jMorris fe Podsiadlowskil 120071 ) or to equatorial wind compression (j Collins et al.lll999l ). Hy- 
drodynamic simulations of the subsequent interaction between the stellar outflow can repro- 
duce th e overall bipolar CSM structure and can possibly explain the origin of the triple-rin g 



nebula (IBlondin fe Lunddqvistlll993l : iMartin fe Arnettlll995 



Morris fe Podsiadlowskil 120071 ). 



At radio fre quencies, an initia l outburst from SN 1987A peaked on day 4 and then 
decayed rapidly (ITurtle et al.lll987h . The radio emission was redetected again in mid-1990 
by the Molonglo Observatory Synthesis Telescope fM OST) (ITurtle et al.l 119901) a nd then 
by the Australia Telescope Compact Array (ATCA) (jStaveley-Smith et al.lll992l ). Long 
term monitoring has shown a monotonic increase in the radio flux since the reappearance 
( jStaveley-Smith et al.l 120071 ). This is believed to be synchrotron radiation from the con- 
tinuous injection of nonthermal particles as the blast wave expands into the CSM. In this 
general picture, electrons are accelerated to ultrarelativistic velocities by the shock and 



- 3- 



then interact with the magneti c field that is amplifi ed by the Rayleigh- Taylor instability 
near the contact discontinuity (jJun &: Norman! Il996l ). Therefore the radio emission is ex- 
pected to be generated between the forward and reverse shocks. A deta iled time-dependen t 
model of diffusive shock acceleration for SN 1987A was first discussed bv lBall fc Kirkl (119921 ). 
Duffy et al.l (119951 ) further improved the model by c onsidering the accelerated ions t hat heat 
the upstream plasma and modify the shock front. I Chevalier fc Dwarkadasl (119951 ) slightly 
modified the physical picture and attributed the radio emission to the interaction between 
the ejecta and an Hn region located inside the optical ring. In this model, the Hn region 
is a shell of swept-up RSG wind photoionized by the BSC The ongoing encounter between 
the shock and CSM also predicts a stea dy increase in X-r ay emission, which has been con- 
firmed by long term X- r ay m onitoring (IPark et al.l 120071 ). Moreover, it was suggested by 
Chevalier fc Dwarkadasl (119951 ) that the Hn region would delay the encounter of the blast 
wave with the optical ring to 2005 ± 3 CE. By then, there was predicted to be a drastic 
increase in the radio flux. 

From an energetic point of view, generation of the radio emission requires much less 
energy compared to other wavebands such as optical and X-ray. Hence, radio data are a 
sensitive tracer of the shock interaction. As the supernova ejecta expand into the CSM, 
the radio morphology will follow the evolution of the forward and reverse shocks and traces 
out the global structure of the dense CSM. High resolution imaging obs ervations with the 
ATCA have been performed since mid-1991 (jStaveley-Smith et al.lll993al ). Fig. [T] shows the 
diffraction-limited images at 9 GHz. With an angula r resolution o f / /9, t hese images already 



reveal t he spatially extended structure of SN 1987A. lBriggsl (119941 ) and lStaveley-Smith et al. 



( 11993a! ) tested and employed the super-resolution to further improve the resolution to C'4 and 
found that the remnant has a limb-brig htened shell mo r pholo gy (Fig. [2]). Applying the same 
technique to observations in 1992-1995. iGaensler et al.l (119971 ) discovered two bright lobes in 
the shell. The lobes were aligned along the major axis of the optical ring and the eastern lobe 
was about 1.2-1.8 times brighter then the western lobe. The ATCA 18 GHz upgrade provided 
twice the angul ar resolution as at 9 GH z, and confirmed all these structures in diffraction- 
limited images ([Manchester et al.ll2005l ). The 18 GHz image is well described by an inclined 
thick equatorial ring with the same orientation as the bipolar CSM. Recently, Tingay et al. 
(2008, in preparation) performed a new VLBI observation of SN 1987A. The 1.4 GHz image 
similarly shows two lobes coincident with the ones in the 9 a nd 18 GHz image s. In the hard 
X-ray band, a thick ring morphology has also been observed (IPark et al.ll2002l ). suggesting a 
similar emission process as at radio frequency. However, t he optical hot spots, which are the 
inwa rd protruding dense gas shocked by the blast wave ( Pun et al. 20021; ISugerman et al. 



20021 ). show no correlation to the radio lobes in position ([Manchester et al.ll2002l ). Although 
the emergence of the optical hot spots coincided with the arrival of the radio-producing shock 
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at those regions (jBall et al.ll200ll ). lManchester et al.l (120021 ) indicate that this had little effect 
on the morphology of the radio remnant. 

As in many astrophysical observation s, quantitative meas urements of the physical pa- 



ramete rs require careful data modeling (e.g. Ng fc Romani|[2008 ). In our case, Staveley-Smith et al 



our case, 

(jl993bl ) demonstrated that a robust way to obtain the remnant size is to model the visibility 
data directly in the u-v plane. They performed two dimensional (2D) Fourier transform of 
a range of models including a spherical shell and a thin ring and t hen fi tted them to the 
visibility data to measure the remnant's expansion. iGaensler et al.l ( 119971 ) employed a thin 
spherical shell model in the fit and found a mea n expansion velocity of ~ 3000 km s~ x be- 
tween 1992-1995. Later. Manchester et al.l (120021 ) improved the modeling by adding two point 
sources to the thin shell model to account for the eastern and western lobes as observed. Re- 



cently, IGaensler et al.l (120071 ) applied the thin shell model to fit all 9 GHz observations from 
1992 to 2006 and suggested a gradual acceleration in the expansion followed by a recent 
deceleration. 

The ATCA imaging campaign has the longest time span for monitoring of SN 1987A at 
any wavelength. With more than 17 years of observations from a single instrument, we are 
now ready to improve the modeling. In this paper, we develop a better model to account 
for more features in the remnant. We will model the observed east- west asymmetry as well 
as the expected light travel time effects. The fits to all the 9 GHz and 18 GHz observations 
from 1992-2008 can then provide a better understanding of the structure and evolution of 
SN 1987A. We describe the ATCA observations and data reduction in §2], then the detailed 
modeling and fitting procedures in $31 The results and discussion are presented in §JH and 
[5j Detailed mathematics of the model is listed Appendix [A] 



2. Observations 

The ATCA is a radio interfer ometer consisting o f six 22-m dishes on a 6-km track, located 



near Narrabri, NSW, Australia dFrater et al.lll992f ). The ATCA observations of SN 1987A 



discussed here were performed at 9 and 18 GHz with 6 km array configurations. At 9 GHz, 
the source was observed approximately every six months with long integration times (6-12 
hr) at two simultaneous frequencies in the range 8.6 to 9.1GHz. In the observatio ns, the 



phase center was placed 10" south of the remnant to avoid any possible artifacts (jEkers 



19991 ). Detailed observational parameters of all the datasets used in this study are listed in 
Table [TJ Two ATCA 18 GHz observations are included in the analysis. Each was taken one 
day before or after a 9 GHz observation, thus allowing a direct comparison between the two 
frequency bands. Except for the last epoch, a preliminary analysis of all the 9 and 18 GHz 
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data has been presented in iGaensler et al.l (120071 ). 



The data reduction was carried out with the MIRIAD package (jSauft et al.lll995l ). We 
examined the data carefully for bad points and rejected scans during poor atmospheric phase 
stability. From 1996 onwards, the source has sufficient signal-to-noise ratio (S/N) to allow 
phase self-calibration, which was performed with a 5 minute solution interval time. This 
time scale is long enough to integrate an appreciable S/N, but is short enough to track the 
atmospheric phase variations. The phase self-calibration results in a significant improvement 
in the fidelity of the complex gains for each antenna. After the calibration, we shifted the 
phase center to the supernova site and then averaged the visibility data in 5 minute intervals. 
Again, this is long enough to obtain good S/N on each baseline, but is short enough to track 
the structure of the source. We note that there is negligible loss in amplitude due to time- 
averaging of the data, because the source is of small angular extent and has been shifted to 
the phase center (IBridle fc Schwab! Il999l ). 



3. Modeling 



Figures [T] and [2] show the diffraction-limited and super-resolved images, respectively, 
of the 9 GHz observations; some epochs were averaged togeth e r to b oost the S/N (see Ta- 
bled]). Following the procedures described by IGaensler et al.l (119971 ). the dirty maps were 
first generated from the calibrated visibility d ata, then deconvolved using a maximum en- 
tropy (MEM) algorithm (jGull fc Danieillll978l ). The MEM model was then convolved with 
a diffraction-limited beam and a super-resolved circular beam with FWHM Pf!4 for Fig. [T] 
and [2J respectively. Finally the models were combined with the unsealed residuals and then 
regridded at a pixel scale of O'.'Ol. The structure observed in the supe r-resolved 9 GHz images 



is confirmed by the 18 GHz diffraction-limited images (see Fig. 1 in Manchester et al.l 12005 



and Fig. [9] of this paper), showing a morphology that is well-d escribed by an inc li ned t hick 
ring. Based on the surface brightness at the center of the ring, Manchester et al.l (120051 ) ar- 
gued that the radio morphology sits somewhere between a spherical shell and an equatorial 
ring. 

Motivated by all these results, we below present a generic inclined torus model to capture 
the overall structure of the radio emission. We emphasize that our model is a simplified 
picture of the physical situation and that we have not attempted to reproduce every single 
detail of the emission. Instead, this study aims to obtain characteristic scales of the emission, 
such as the size and height of the emitting region, and the evolution of these parameters 
with time. 
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3.1. Equatorial Belt Torus Model 

The extra-planar structure of SN 1987A mentioned above obviously does not fit into the 
simple picture of a spherical shell or a 2D ring, but is better described by a 3D torus with a 
finite thickness. To incorporate these features, we propose an equatorial belt torus model, i.e. 
a spherical shell truncated at arbitrary latitude. This gives a continuous variation between a 
ring and a spherical shell geometry. Fig. [3] illustrates this model. It is a truncated optically 
thin shell with uniform emissivity and inclined to the observer line of sight. The geometry is 
specified by the radius R (averaged between the inner and outer shells), the shock thickness 
5 (as a fraction of R) and the half opening angle 9. Additional parameters include the flux 
density f , center position (x ,y ), inclination ( and position angle (P. A.) ^ on the plane 
of the sky. In the limiting case of 9 = 0°, the torus is reduced to an equatorial ring, while 
9 = 90° gives a spherical shell. The east-west asymmetry of the remnant is modeled by a 
linear gradient in the torus emissivity across the equatorial plane (not in the plane of the 
sky). Hence two more parameters are introduced: the magnitude g and the P. A. <fi of the 
gradient. The former is normalized to the fractional change in flux density at R, and the 
gradient direction is defined from the brightest point towards the faintest point, with the 
P. A. measured from the projected observer line of sight (see Fig. [3]). Finally, we correct for 
light travel time effects. Since the remnant flux is increasing, the north rim, which is about 
a light year closer to the Earth than the south, appears brighter. Also due to the finite 
speed of light, expansion of the remnant distorts the torus and shifts its center to the north. 
To derive corrections for both effects, we employed linear approximations to the expansion 



and light curve 



ocally at e ach epoch, with the rate of c hange determined from the results of 



Gaensler et al.l (120071 ) and IStaveley-Smith et al.l (120071 ). We found that the north rim is at 



most 20% brighter than the south as a result of the light travel effects. On the other hand, 
the effects of expansion are only of second order, at the level of less than 1%. Therefore, the 
latter will be ignored throughout the following analysis. A detailed mathematical description 
of the torus model and light travel travel time effects is given in Appendix |A] 



3.2. Fitting Procedure 



Staveley-Smith et al.l (Il993bl ) showed that a robust way to estimate the remnant's ge- 



ometry is to model the visibility data directly in the u-v domain, for which the measurement 
errors are uncorrelated and the fitting is not complicated by the deconvolution process. In 
order to do that, we first have to derive a 2D Fourier transform for the inclined torus model. 
This was achieved by building up the model torus from t hin rings, which can be transformed 
to Fourier space analytically via the Hankel transform (IBracewelll 120001 ) . The linear gradi- 
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ents in the model, including the asymmetry and light travel time correction, are transformed 
using the convolution theorem and are turned into derivatives in the u-v domain. Putting 
these together, we obtain an analytic expression of the torus model under a Fourier transform 
(see Appendix [A] for the detailed formulae). 

Once the model is established, the fitting procedure is by the minimization of 



X 



E 



model (w, v) 



a(u, v) 



where a is the statistical uncertainty associated with each measurement. Given that the rms 
noise levels in our final images is close to the thermal noise of the antenna receivers, we argue 
that other sources of errors in the measurement, including the systematic errors due to the 
calibration or instrumental effects are negligible. The \ 2 minimization is carried out by the 
task uvf it in MIRIAD, with modifications to the code for the torus model. As discussed, 
the symmetry axis of the system ha s been well measured from optica l observations of the 
triple-ring system and of light echoes (jSugerman et al.ll2005bl ; |Punll2007l ). Therefore, we fixed 
the o rient a tion o f the torus during the fit, by adopting the values \& = — 7?6 and ( = 43? 4 
from iPunl (120071). The uncertain ties in the fits were estimated through a bootstrapping 
method (jEfron fc Tibshiranilll993l ). since the errors in the fits are likely to be non-Gaussian. 
Following the procedure described in lPress et al.l (119921 ). boostrap samples were generated by 
randomly drawing the visibility measurements from an observation, with replacement. After 
the same number of visibilities as the data have been drawn, the bootstrap replica is fitted 
with the same procedure to obtain the torus parameters. For each observation, we fitted 
2000 bootstrap iterations, and the 68% intervals for each parameter's distribution are quoted 
throughout as the la uncertainty. This gives the projection of multi-dimensional errors that 
accounts for the correlations between fitting parameters. As the S/N in the observations 
increases towards the later epochs, the parameter uncertainties due to measurement errors 
become negligible and the systematic errors in the modeling dominate. In the real world, the 
detailed structure of SN 1987A is much more complicated than what our simple model can 
capture, resulting in systematic errors in the fits. We have attempted to quantify this term 
for the radius, which is the parameter we are most interested in. Since the exact structure 
of the remnant is unknown , the radius measurement is always model-dependent. Similar 
to Manchester et al.l (120021 ). we fitted the 9 GHz observations to various plausible models, 



including a spherical shell and a shell plus two point sources. In general, these fits obtain 
slightly different values of R and hence the variation provides a handle of the systematic 
errors in the modeling. 
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Results 



4.1. 9 GHz Observations 



The best-fit torus parameters are listed in Table [2] and plotted in Fig. Most of the 
fits to the ATCA 9 GHz data are robust, except for a few epochs with very large \ 2 values 
that are possibly due to poor calibration. The fits are generally better towards the late 
epochs, however as shown in Table [2j the typical reduced \ 2 value is ~ 1.5. This is still 
relatively large, suggesting that our model is not a perfect fit to the observations, and that 
there are unmodeled features in the data which contribute to the systematic residuals. This 
result is not unexpected since the model aims to capture only the characteristic features of 
the emission, rather than reproducing every small-scale structure. Nevertheless, our model 
does provide a bette r description to the data than a simple thin spheric al shell model used 
in previous studies (jGaensler et al.l 119971 . 120071 : Manchester et al.ll2002l ). Fig. @] plots the 
visibilities of the data and models along two orthogonal axes in the u-v domain. The differ- 
ence between a torus and a spherical shell models is obvious, in particular along the w-axis, 
because this direction corresponds to east-west in the sky, where the radio emission shows 
large asymmetry. The plots clearly indicate that the asymmetric torus model is a better fit 
over a simple shell, and, at least to the first order, well describes the east-west asymmetry. 
On the other hand, w e note that a similar plot with azimuthally averaged u-v data (Fig. 6 in 
Gaensler et al.l 120071 ) masks the discrepancy between the data and a simple shell model. To 
provide a quantitative comparison between different models, we have repro duced the fits with 



a thin shell model and a thin spherical shell plus two points model as in Manchester et al. 



( 120021 ). Table [3] shows that the torus fit always provides lower x 2 / u values than the shell 
fit. The shell plus two points model and the torus model have comparable goodness-of-fit, 
but the former requires two more fitting parameters. Therefore, we conclude that our torus 
model is a substantial improvement over the simple models used previously. 

Although we did not attempt to provide a physical model to interpret the observed 
radio light c urve, we note that the fluxes obtained from the fits are consistent with the values 
reported by IStaveley-Smith et al.l (120071 ). In the following discussions, we will concentrate 
on the torus geometry. 

Fig. [6] shows the best-fit model for epoch 2008.0, with the left panel illustrating the 
model at very high spatial resolution. The images were generated in the u-v domain with 
complete uv coverage, then Fourier transformed back to the image plane. The right panel 
shows the image convolved with a Of! 4 FWHM Gaussian, the same as the super-resolved 
beam used in Fig. [2J The inclined torus model has successfully captured the aspect ratio of 
the overall structure, and has also reproduced the two lobe features, thus confirming that 
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they are due to projection of the optically thin equatorial belt. 

A better way to compare between the model and data is to generate an image of the 
model with identical u-v sampling as the corresponding observation. Then the sampling 
effects will appear identically in both images, allowing a direct comparison. Fig. [7] & [8] show 
the images of the data compared to the best-fit models for selected epochs. The models were 
generated with the approach described above and then imaged using the same procedures 
as mentioned in §|2J However, since the MEM algorithm is not applicable to images with 
negative values, we did not attempt to deconvolve the residual maps. We emphasize that the 
residual maps presented here only illustrate a qualitative comparison between the data and 
the model, rather than a detailed quantitative description of the underlying difference. The 
most obvious feature in the residuals is the flux deficit of the model torus in the southeast and 
northwest, and the excess emission at the center. The discrepancies are due to limitations 
of our relatively simple model. As indicated by the super-resolved images, the eastern lobe 
in the data is not reflection symmetric about its brightest point, but extends further to the 
south. This cannot be easily described by the model, which approximates the asymmetry by 
a linear gradient. Although less obvious, a similar case occurs for the western lobe as well. 
On the other hand, the excess flux at the center of the model torus is due to the abrupt 
cutoff of the high latitude emission in the model. Since the inclination angle ( is fixed during 
the fit, the ellipticity depends solely on the opening angle 9. A large 9 is required by the fit 
in order to match the overall geometry, hence resulting in excess emission filling the interior. 
In real life, the radio emissivity is likely to be latitude dependent and to decay smoothly at 
high latitude, providing a slightly different overall ellipticity. A more realistic model should 
incorporate all these effects to address the above discrepancies. 



4.2. 18 GHz Observations 

The best-fit model for the 18 GHz observation of 2003 July is shown in Fig. [91 with the 
data and model visibilities plotted in Fig. [TDJ The diffracted-limited images were deconvolved 
the CLEAN algorithm. Similar to Fig. [7] & El no deconvolution has been applied to the 
residual map in order to avoid running the algorithm into noise. As shown in Table |2j the best 
fits have small x 2 l v ~ 1? which are due to large measurement uncertainties (see Fig. fTOl) . The 
best-fit parameters are plotted by the open circles in Fig. They are all consistent with the 
9 GHz results, except that 9 is slightly smaller at 18 GHz. While the discrepancy in 9 seems 
to suggest a frequency-dependent emission height that implies different physical conditions 
between the equatorial plane and at higher latitude, it could also be the result of systematic 
errors in the fits which we do not fully understand. As the remnant is getting brighter, future 
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observations will help to investigate this problem. Following iManchester et al.1 (120051 ). slices 
through the 18 GHz diffraction-limited images along different P.A. were plotted in Fig. [TTJ 
In general, the two peaks in the model match well to the observation. The only exception 
is at 150°, where the data are ~ 25% brighter than the model in the east. As discussed in 
§4. II above, this is attributed to the unmodeled southern extension of the lobe. The plot also 
illustrates the excess emission near the center of the model torus, again consistent with the 
9 GHz results. 



5. Discussion 



5.1. Remnant Size 



In general, the size of the to rus is well-determined from the fits and is consistent with that 



deter mined in previous studies (IGaensler et al.lll997l ; IManchester et al.ll2002l ; iGaensler et al. 



20071 ). The best-fit radii from the torus and shell fits are comparable, but those from the shell 
pl us two points mod e l are always slightly larger. A similar discrepancy was also reported 
by IManchester et al.l (120021 ). suggesting that the radius measurement is somewhat model- 
dependent. The average discrepancy is ~ 5% for the later epochs, which is attributed to 
the systematic errors in the modeling as mentioned in §3.21 In the following discussion, we 
assigned a 5% systematic errors for R. Hence the quadrature sum of these and the statistical 
errors obtained from bootstrapping provides the total uncertainties in R, which are shown 
by the dotted error flags in Fig. 



Compared to the measurements published for the X-ray shell (jPark et al.l 120071 ) , the 
radio torus is larger by 1 0-15% which is too la rge to be accounted for by the measurement 
uncertainties. However, IGaensler et al.l (120071 ) analysed the Chandra data in the Fourier 
domain and found that the X-ray and radio shells agree to better than 1% in size. Hence 
they argued that the discrepancy is due to different measuring t echniques, rather than any 
physical difference. At optical wavelengths, iMichael et al.l ( 119981 ) reported that the reverse 
shock was located at 70% of the inner circumstellar ring around day 3900. This gives a 
radius of 0'.'64, which is slightly s maller than the torus size at that epoch. From the steady 
increase in the radio light curve, IManchester et al.l (120051 ) argued that the radio-emitting 
electrons could be filling the region between the shocks. Although the thickness 5 is not 
too sensitive in the fits and generally very small, this is still consisten t with the prediction. 
With a density profile p{r,t) oc t~ 3 (r/t)~ 9 in the supernova ejecta (jEastman fc Kirshner 
19891 ). numerical simulations suggest that during the early stage of the SNR evolution, the 
forwa rd and reverse s hocks are separated only by 10-20% in radius (e.g.. iTruelove fc McKee 
1999h . For SN 1993 J, iBartel et al.1 J2007h reported a shell thickness of 30% of the inner shell 
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radius. Despite the above results assuming a spherical symmetry, if 5 in our case is of a 
similar order, given that the whole remnant is barely resolved i n the 9 GHz observatio ns, 5 
likely remains unresolved. This also supports the conclusion of iGaensler et al.l (120071 ) that 
both thick and thin shell models provide equally good fits to the data. 



5.2. Rate of Expansion 



At the earliest epoch, the remnant has a radius of about 0'.'6. As lGaensler et al.l (119971 ) 
pointed out, this implies a minimum mean expansion rate of 30000 km s _1 during the period 
1987 to 1992, in the free expansion phase. This expansion was then followed by a drastic de- 
celeration to ~ 3000kms _1 at around the time the radio emission re-emerged (IGaensler et al. 



19971 ). To determine the expansion since day 1800, we employed a linear fit to R for all the 



9 GHz observations and obtained an expansion velocity of v — 4000 ± 400 km s -1 over 
the period 1992-2008. The small x 2 / u = 0.2 in the fit indicates that the measurement 
uncertainties are over-estimated. On the other hand, if only the statistical errors are con- 
sidered, the result is similar (v = 4150 ± 30 km s _1 ), but with smaller errors and a much 
larger x 2 jv = 9.0. IGaensler et al.l (120071 ) performed a quadratic fit to the remnant size and 
suggested an acceleration in the expansion. We also found a similar trend in our results. 
However, a F-test suggests that the quadratic fit is not statistically significant: there is 
a 20% probability that the smaller \ 2 i n the quadratic fit would occur solely by chance. 
Therefore, we conclude that the simpler model (i.e. linear fit) is preferred and adopt the 
linear expansion v = 4000 ± 400 km s _1 in the following disc ussion. This va l ue lie s be- 
tween v = 4700 ± 100 km s _1 and ~ 3700 km s -1 obtained by IGaensler et al.l (120071 ) and 



Manchester et al.l (120021 ) respectively using a thin shell model, and agree s with the prediction 
v = 4100 km s _1 from hydrodynamic simulations (IBorkowski et al.lll997l ). It is also similar to 
the ex pansion rate of the reverse shock at v = 3700 ± 900 km s _1 reported by iMichael et al. 



( 119981 ). On the other hand, the X-ray expansion is apparently much faster before day 6200. 
Early measureme nts showed an X-ray expansion of v = 4200 km s _1 (IPark et al.l 120041 ) . but 
Park et al.l (120071 ) recently suggested acceleration to a higher velocity of v ~ 6000 km s" 1 
over the period covering 5000-6000, before the expansion decelerates to v — 1400 km s _1 
near day 6200. This could indicate that the X-ray emission is located differently as the radio 
emission and much closer to the forward shock. 

A comparison between the radio and optical observations suggests that the supernova 
blast wave may have started to encounter the dense CSM in the past few years. The op- 
tical ring, the emission from which represents fluorescence after ionization of the CSM by 
the supernova UV flash, indicates the inner boundary of the dense swept-up CSM. As seen 
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by the HST, th e central equatorial ring has a radius of 0'/81-0'/85 (IBurrows et al.l 11995 



Plait et al.lll995l ). Accounting for the uncertainties in both the radio and optical measure- 



ments, our results suggest that the radio emission could have reached the optical inner ring 
around day 6000-6500, hence the forward shock could have started t o expand into the dense 
CSM a s early as 2004. This falls nicely within the predict ions by IChevalier fc Dwarkadas 
(119951 ). iGaensler et ajj (|1997l) and Manchester et al. d2002h . but disagrees with other esti- 
mates such as those of lLuo et al.l (11994 ) and lBorkowski et al.l (119971 ). The X-ray observations 



also provid e evidence 



6000-6200 ( Park et al. 



or th e encounter. The rapid increase in the soft X-ray flux near day 



20051 ) and the drastic deceleration of the X-ray expansion around day 



6200 (IPark et al.l 120071 ) reveal that the forward shock has encountered a steep gradient in 
the ambient density. We expect a similar deceleration in the radio expansion to occur very 
soon. 



5.3. Emission Height 



In addition to the radius, the opening angle 9 is also an important parameter that probes 
the CSM geometry. The best-fit values of 9 shown in Fig. [5] are substantial, indicating a 
significant level of radio emission from mid-latitudes. This extra-planar structure implies the 
spread of dense CSM over large regions on either side of the equatorial plane, and suggests 
that the prominent optical ring is merely the densest part of a more complicated 3D structure. 
In the later epochs, 9 stays nearly constant around 40°. This gives an inferred emission height 
up to ~ 5 x 10 17 cm from the equatorial plane. Our result ind i cates a radio emission geometry 
very similar to that of the reverse shock. Michael et al.l (119981 . 120031 ) reported a rapid decrease 
of Lya and Ha emission beyond 40° and hence deduced that the reverse shock is mostly 
confined within 30° of the equatorial plane. Our fitting results also provide the characteristic 
structure of the CSM, which c an be compared to hydro dyn amic simulations such as those o f 
Blondin fc Lunddqvistl Jl993h . Martin fc Arnettl Jl995h and Morris fc Podsiadlowski (l2007h . 
In those simulations, the exact shape of the contact discontinuity between the BSG and RSG 
winds is highly model-dependent. So me models suggest a rath er constant opening angle over 
a large range of radius (e.g. Fig. 3 of Martin fc Arnettl Il995l ). which is similar to our case. 
Nevertheless, all the simulations indicate that the overall geometry of the swept-up CSM is 
peanut-like, with a dense equatorial ring at the waist. This picture predicts that the radio 
emission will extend with time to higher latitude. Although the expected trend is not found 
in the current data, we anticipate 9 to increase significantly in future observations according 
to the above models. Eventually the shock will reach the dense CSM in the polar region and 
the radio morp hology will be more c ircular. Our results do not show any polar protrusions 
as suggested by lBlondin et al.l (119961 ). If the forward shock continues to expand at a similar 
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rate, it may take another 10-20 years to reach the polar region before we can verify these 
predictions. As a note, Fig. [5] suggests some hints that 9 is slightly larger in the early epochs. 
If this is real, it may indicate the encounter of ejecta with dense regions at higher latitude 
at the early stages. Along with the shock expansion, the radio emission is then gradually 
dominated by other regions closer to the equatorial plane, resulting in a gradual decrease in 
6. 



5.4. Asymmetry 



The super-resolved images in Fig. [2] clearly show that the eastern lobe is brighter than 
i ts western counterpart. This feature has also been rep orted in the X-ray observations 
(IBurrows et al.ll2000l ; iPark et al.ll2002l ; iMichael et al.ll2002l ). In our torus model, the asym- 
metry is quantified in terms of a linear gradient with direction and magnitude g. Fig. [5] 
shows a trend of increase in prior to 1998. Although this occurred at a similar time as the 
optical hot spots first appeared, comparing to the hot spot positions shows tha t they do 
not coincide. Therefore we arrive the same conclusion as Manchester et al.l (120021 ) that the 
detailed morphology of radio and optical emissions are not correlated. The only similarity is 
that they both brightened first in the eastern half, indicating a global asymmetry along the 
east-west direction. On the other hand, the variation in may be associated with the slight 
decrease in 6 around that period. This could suggest that the blast wave has encountered 
different dense regions before and after ~ 1998. In the late epochs, stays around 110°, 
which corresponds to a P. A. of 96° on the plane of the sky. This m atches the bright X-ray 
lobes located at P.A. 90° and 270° as reported bv IPark et all (|2002h . 



Although there were indications from earlier data that the asymmetry increases with 
time (IGaensler et al.lll997l ; Manchester et al.ll2002l ). the longer data span presented here does 
not support this. From 1995 onward, g stays at a constant level of ~ 40%. This gives a flux 
ratio of 1.5 between the two halves of the torus at the equatorial plane. The projection onto 
the plane of the sky reduces this ratio slightly. Thus as seen by an observer, the eastern 
lobe appears to be ~ 1.4 times brighter than the west. In order to confirm this result, we 
compared to the degree of asymmetry inferred from the shell plus two points fits. The surface 
brightness ratio between the eastern and western lobes is listed in Table [3] again this does 
not show any significant long-term variations, but suggests a rather constant flux ratio of 
about 1.3, which is consistent with the torus estimate. Hence we conclude that the surface 
brightness of the eastern half is consistently 30-40% larger than the west. This is very close 
to the 40% asymmetry in the reverse sho ck (IMichael et al.ll2003T) . also similar to the 15-50% 
asymmetry in the 0.3-6 keV X-ray band (IMichael et al.ll2002l ). While the exact origin of the 
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east-west asymmetry still remains unclear, it seems likely that the asymmetry either lies 
in the initial velocit y of the ejecta, which would imply an asym metric explosion, or in the 
density of the CSM (IGaensler et al.lll997l ; ISugerman et al.ll2002l ). In principle the evolution 
of the radio morphology can help distinguish these two scenarios. If the eastern lobe expands 
faster, it will cause a systematic shift of the torus center towards the east. Unfortunately, 
the phase self-calibration process removes the absolute astrometric information in the data, 
precluding a direct comparison between epochs. Nonetheless, the shell plus two points fit 
can provide some insights into the relative positions of the features. Although this model 
is less preferable than the torus model, both of them obtain a comparable goodness-of-fit, 
therefore we believe that the former can still provide a sensible description of the location of 
the lobe centroids. As listed in Table [3X the offsets of the point sources from the shell center 
increase at a similar rate as the shell radius. This would not be the case if the CSM were 
denser towards the east, since the higher density would lead to a larger deceleration for the 
eastern lobe. In the other scenario, if the ejecta in the east had a higher velocity, then it 
would arrive at the equatorial ring first and would start to decelerate. Hence we expect the 
ratio of offsets to vary in future observations. However, we emphasize that this is only a very 
rough estimate. In particular, if the shell plus two points model is not a good description of 
the physical structure, excess flux from the lobes will shift the center of the underlying shell, 
introducing large systematic errors in the positions. 



5.5. Limit on Radio Flux of a Central Source 



In spite of extensive searches in radio and optical wavelengths, no pulsar has yet been 
found within the supernova remnant ( Manchesterll2007l ) . Our fitting scheme can also provide 



a detection limit on any possible central object. We applied the same fitting procedure using 
the torus model with a point source fixed at the center. The flux density of the point 
source was held fixed during the fit. We adjusted its level until the change in the best- 
fit x 2 value exceeded a certain limit. To begin with, we applied the fit to the 1996 July 
observation, in which the remnant has high S/N but is still not too bright to preclude a 
sensitive search. The fit gives a 3a upper limit of 0.3 mJy Applying the fit on several other 
epochs obtained very similar results. Hence we conclude that this is the 3a detection limit 
on any unresolved central source, including a neutron star or a compact pulsar wind nebula, 
for the ATCA observations at 9 GHz. As a comparison, th e most sensitive li mit on pulsation 
at this frequency yields an upper flux limit of 0.058 mJy ([Manchester! 120071 ) . 
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5.6. Comparison to Other Radio Supernova Remnants 



Among about ten resolved radio supernovae (RSNe) (see iBietenholzl 120051 ). SN 1987A 
shows a very different characteristics in the asymmetry, morphology and expansion as com- 
pared to the others. For example, SN 1993J in M81, which is the best studied RSN besides 
SN 1987A, was found t o be circular with i n 3% over 7 years of VLBI observations. Based on 
a statistical argument, iBietenholz et al.l (120031 ) suggested that the emission is intrinsically 
spherical. On the other hand, only two other RSNe show large deviation from a circular shell. 
VLBI imaging of SN 1 9 86J in NGC 891 revealed a distorted shell structure with a hot spot 
(IBietenholz et al.ll2004j ). iPerez- Torres et al.l (120021) suggested this could be d ue to protrusions 



in th e progenitor's wind (IBlondin et al. 



19961 ). or jet-induced core collapse (jKhokhlov et al. 



19991). Another example of a non-spherical RSN is SN 2004et in NGC 6946. A recent study 
by iMarti-Vidal et al.l (120071 ) interpreted the asymmetry as either a shell with two hot spots, 
which could be similar to the case of SN 1987A, or an anisotropically expanding shell with 
a protrusion. 

The main distinction between SN 1987A and other RSNe lies in its rate of expansion. 
As discussed, the expansion of SN 1987A must have undergone a rapid deceleration in the 
first 1800 days. In contrast, for all other RSNe with detected expansion, their sizes are 
well-approx i mated by a power-law R oc t m , with the deceleration parameter m ~ 0.7 — 1 
(IBietenholzl 120051 ) . The most extreme example is SN 1979C in M100, wh ich has shown 
nearly free expansion (m ~ 1) for over 22 years (IBartel &: Bietenhola 120031 ) . The unusual 
morphology and evolution of SN 1987A can be largely attributed to its complex circumstellar 
environment. Thus, this extraordinary event provides a unique example in the study of 
supernova shell dynamics, stellar wind interactions and stellar evolution. 



6. Conclusions 

In this study, we have developed a scheme to model the interferometeric u-v data for any 
3D object with axisymmetric structure. An application to ATCA 9 and 18 GHz observations 
of SN 1987A has improved the spatial modeling of the radio remnant. An inclined equatorial 
belt torus model in 3D was employed, with the incorporation of light travel time effects and 
of an overall east-west asymmetry. The projected brightness distribution of the model is 
transformed into 2D Fourier space analytically, then is fit directly to the visibility data. 
The torus model successfully captures the characteristic features of the radio morphology 
and provides substantial improvements over the simple spherical shell fit used in previous 
studies. The best-fit torus geometry shows significant extra-planar emission extending to 
mid-latitudes ~ 40°, and a 40% flux difference between the eastern and western lobes. 
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Evolution of the shell radius with time implies a constant expansion rate of 4000 ±400 km s" 1 
over the period 1992-2008 covered by the observations; the radio remnant is probably now 
starting to encounter the optical inner ring. Our results suggest a torus geometry and 
expansion speed both very similar to those seen in Ha and Lya emission, indicating that the 
radio synchrotron emission originates between the supernova forward and reverse shocks. 
In summary, the torus fitting provides a powerful tool to probe the 3D structure of the 
radio emission, showing that the overall picture of the dense CSM is consistent with the 
hydrodynamic simulation of the progenitor star's interacting winds. Future evolution of 
the radio morphology can help differentiate between specific evolutionary histories for this 
system, and can address the apparent east-west asymmetry in the ejecta distribution. 
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A. Two Dimensional Fourier Transform of an Inclined Torus Model 



In this section, we give a mathematical description of an inclined torus model and 
derive an analytic expression for its Fourier transform. Although a specific example of a 
torus model is considered here, the expression could be easily generalized to any objects 
with axisymmetry. 

We first set up a coordinate system (x', y', z') with z' along the symmetric axis and y' 
along the projected line of sight at the equatorial plane. In this frame, a model torus (or 
any axisymmetric object) can be built up from rings centered at (0, z') with radius r 



I torus (x',y') 



w Ii^f dr dz' 



where the intensity of each ring is given by J"°, g = fS^r— a/ x' 2 + y' 2 J . In our case, the 
weight fact or has value w = 1 as the torus has a uniform emissivity. Applying the Hankel 
Transform (IBracewelll [2000J ) . each ring in the Fourier domain is given by 



wing. 



V? + V 2 



where Jo is the Bessel Function of the First Kind. In order to account for the asymmetry 
observed in the remnant, we assume that the flux density of the torus varies linearly at a 
rate g along the direction <fi in the equatorial plane (see Fig. [3]). Hence the intensity of each 
ring is modified to 

jasym = ^ ^ ^ + y , ^ ^ + ^jring _ 



With the Convolution Theorem (IBracewelll 120001 ) . the Fourier transform of an asymmetric 
ring is then 

F asym (u,w) = ^{/ asym } 

= T{g(x' cos (j) + y' sin 0) + 1} * *F{I Ting } 

ig_ 

2vr 

ig 



cos 4>8'{u)5{v) + sin <j)5{u)5'{v)} + 5{u)5{v) } * *F rms {u,v) 



cos 0— + sm <p— | + 1 
Z7T \ ou ov 



{U,V} 



irg 



asm 



<p + v cos 4>)J\ (2ix rVu 2 + v 2 j + Jo ^27r rVu 2 + v 2 j , 



This preprint was prepared with the AAS IATgX macros v5.2. 
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where ** denotes a 2D convolution such that 



F(u, v)G(u — u',v — v') dv! dv' 



for any 2D functions F(u,v) and G(u,v). 5(u),S(v) are Dirac's 5 functions and = j- 
and ) = 



dv ' 



The model torus has to be projected onto the plane of the sky before it is compared to the 
observations. Here we set up another coordinate system (x,y, z), with z along the observer 
line of sight and y along the projected torus axis on the image plane. The transformation 
between the two frames is given by 



x 

y 



X 



y' sin ( — z' cos £ 
y' cos ( + z' sin ( 



where ( is the inclination angle between the torus equatorial plane and the observer line of 
sight (see Fig. [3]). The projection essentially compresses the rings into ellipses in the image 
plane, with center at (x, y) = (0, —z' cos and semi-minor axis r sin (. Hence an asymmetric 
ring projected on the image plane can be expressed as 



g Ixcos0+ ^sin°c C sin ' 



+ 1 



fSlr 



x- 



y+z' cos £ 
sin £ 



Following the scaling and shifting properties of the Fourier Transform (IBracewelll 120001 ) . for 
any constants a and b, 



F{f{x -a,y- b)} 
F{f{ax,by)} 

hence 

F pT °i(u } v) = T{P Toi } 
up to an overall normalization factor. 



-i2n(au+bv) 



; iW ^ cosC F asym (M,t;sinC) 



To correct for light travel time effects, we assume the torus flux density / increases 
linearly during the light travel time across the remnant, i.e. / m ft + /o, for constants / and 
/o- While the latter is a fitti ng parameter, the former is determined from the local derivative 
of the light curve given by IStaveley-Smith et al.l (120071 ). The relative light traveling time 
from each point of the ring to the observer is 



a + z d ( > a , ' • t\ d ( V 

At = - = - (y cos ; + z sin Q = - : j 

c c c \tan(. 



+ 



sinC 
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where d — 51 kpc is the distance to LMC, and z is in the units of radians. Hence the flux 
density at each point is 

df ( y z' 



f 



c \tan£ sin( 
and the torus intensity is modified by a linear gradient of 



+ /o 



df f y + z> 



c f \tan£ sin£ 



+ 1 



Employing the Convolution Theorem again, we obtain the Fourier transform of an asym- 
metric ring including a light travel time correction 



F light {u,v) 



1 d 



df_ 

cf ^2-71" tan £ <% ' sin( 



F pr °i(u,v) . 



Finally, the model is shifted and rotated according to (xo,yo) and \1> respectively, 
F final (w, v) = e -iMxov+vov) F u&t( u cos # _ v S i n ^ u s i n y + v cos 

and integrated through r and z' since the Fourier transform is linear. 



F tOTUS {u,v) 



wF^f dr dz' 
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Fig. 1.— Diffraction-limited 9 GHz ATCA images of SN 1987A over the period 1992 to 2008. 
Some epochs have been averaged together to boost the S/N (see Tabled]). The scale is linear, 
ranging from -0.2 to 31mJybeam _1 . The contours are at levels from 5 to 30mJybeam _1 
with 5mJybeam _1 intervals. The synthesized beam for the first epoch is shown in the first 
panel as an ellipse. The beams for the other epochs vary slightly in orientation, but are 
similar in size. 
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Fig. 2.— Super-resolved 9 GHz ATCA images of SN 1987A over the period 1992 to 2008. 
Similar to Fig. [U some epochs have been averaged to boost the S/N (see Table [1]). The 
scale is linear, ranging from -0.2 to 10.8 mJy beam -1 . The contours are at levels from 2 to 
lOmJybeam -1 with 2mJybeam _1 intervals. The synthesized beam, which is shown in the 
first panel, is a circular Gaussian with FWHM / .'4. 
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Fig. 3. — Schematic diagrams illustrating the geometric parameters in the torus model. Left: 
cross-section of the torus model, as viewed perpendicular to the observer line of sight. Middle: 
the torus in the equatorial plane showing a linear gradient, defined from the brighter side 
toward the fainter side. Right: A model torus, which is generated with parameters similar 
to that of the best-fit model for the most recent epoch (2008.0), projected onto the image 
plane (i.e. the plane of the sky). 
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Fig. 4. — Amplitude of the visibility data in the most recent epoch (2008.0) of ATCA 9 GHz 
data compared to the best-fit torus and thin spherical shell models. The torus and shell 
models are plotted by the solid and dotted lines respectively. The two panels show slices 
through the u-v domain along the axes u — and v — 0. 
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Fig. 5. — The best-fit flux density (/), radius (R), opening angle (9), shock thickness (5), 
gradient angle (0) and gradient magnitude (g) at each epoch, obtained from the equatorial 
belt torus model as described in §3.11 The results for the 9 and 18 GHz observations are 
plotted as filled and open circles respectively. Note that we did not plot the flux density 
for the 18 GHz observations, since the absolute flux scale has not been fully calibrated. The 
uncertainties in flux density is negligible. The dotted error flags in R indicate the total 
uncertainties, which are the quadrature sum of the statistical and systematic errors. The 
best-fit linear expansion in radius of 4000kms _1 is shown by the straight line. 




Fig. 6. — Left: The best-fit model for the most recent epoch (2008.0), with the parameters 
listed in Table O Right: The same model convolved with a 0'.'4 FWHM Gaussian kernel, as 
the super-resolved beam used in Fig. [2] The brightness scales in these images are arbitrary. 
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Fig. 7. — Super-resolved images at 9 GHz (left column) compared to the best-fit models 
(middle column) for selected epochs. These images were generated with the same procedure 
described in £j3l The synthesized beam, which is a circular Gaussian with FWHM 0'/4, is 
shown in the first panel. The right column shows the dirty maps of the residual visbili- 
ties, for which no deconvolution has been applied. The scales are linear, ranging from 
to lOmJybeam -1 in the left and middle panels and from -1.5 to 1.5 mJy beam -1 for the 
residuals. 
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Fig. 8. — Super-resolved images at 9 GHz (left column) compared to the best-fit mod- 
els (middle column) for the last epoch (2008.0). The scales are all linear, ranging from 
10 to 40mJybeam _1 for the diffraction-limited data and model images, and from to 
lOmJybeam -1 for the super-resolved images. As in Fig. [71 the dirty map of the residu- 
als (right column), which are not deconvolved, has a linear scale ranging from from -1.5 to 
1.5 mJy beam" 1 . The synthesized beams are shown by the open ellipse and circle in the left 
column. Note that the model images shown here have been generated with identical u-v 
sampling as the data, while Fig. [6] was obtained by direct Fourier transform of the model 
visibility assuming a complete u-v coverage. 



-30- 



DaLo. 



Model 



Residual 




j i i i i i i i i_ 



1 1 

RA offset (arc see) 



1 -1 

RA offset (arcsee) 



1 -1 

RA offset (aresec) 



Fig. 9. — 18 G Hz diffraction-limited i mages of the data, model and the residual map for 
epoch 2003.6 ([Manchester et al.ll2005l ). The synthesized beam is shown in the left panel. 
The linear scales range from to 3.5 mJy beam -1 in left two panels and and from —0.5 to 
0.5 mJy beam -1 for the residual map. Note that no deconvolution has been applied to the 
latter. 




Fig. 10. — Similar plot as Fig. H] for the 18 GHz observation at epoch 2003.6. Note that the 
left panel is actually a slice at 10° to the f-axis, due to the u-v coverage of the observation. 
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Fig. 11. — Slices through the 18 GHz diffraction-limited images in Fig. [9] at different position 
angles. The offset is radial from the torus center and is positive toward north and/or east. 
The data and model are shown by the solid and dotted lines respectively. The zero levels 
for each slice are shown by the dashed lines. 
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Table 1. Observational parameters for the datasets used in this study. 



Observing date Days since Array Frequency Time on Epoch shown 

supernova Configuration (MHz) source (hr) in Fig. [T] & [2j 



9 GHz Observations 






1992 Jan 14 


1786 


6B 


1992 Mar 20 


1852 


6A 


1992 Oct 21 


2067 


6C 


1993 Jan 4 


2142 


6A 


1993 Jan 5 


2143 


6A 


1993 Jun 24 


2313 


6C 


1993 Jul 1 


2320 


6C 


1993 Oct 15 


2426 


6A 


1994 Feb 16 


2550 


6B 


1994 Jun 27-28 


2681 


6C 


1994 Jul 1 


2685 


6A 


1995 Jul 24 


3073 


6C 


1995 Aug 29 


3109 


6D 


1995 Nov 6 


3178 


6A 


1996 Jul 21 


3436 


6C 


1996 Sep 8 


3485 


6B 


1996 Oct 5 


3512 


6A 


1997 Nov 11 


3914 


6C 


1998 Feb 18 


4013 


6A 


1998 Feb 21 


4016 


6B 


1998 Sep 13 


4220 


6A 


1998 Oct 31 


4268 


6D 


1999 Feb 12 


4372 


6C 


1999 Sep 5 


4577 


6D 


1999 Sep 12 


4584 


6A 


2000 Sep 28 


4966 


6A 


2000 Nov 12 


5011 


6C 


2001 Nov 23 


5387 


6D 


2002 Nov 19 


5748 


6A 


2003 Jan 20 


5810 


6B 


2003 Aug 1 


6003 


6D 


2003 Dec 5 


6129 


6A 


2004 Jan 15 


6170 


6A 


2004 May 7 


6283 


6C 


2005 Mar 25 


6605 


6A 


2005 Jun 21 


6693 


6B 


2006 Mar 28 


6973 


6C 


2006 Jul 18 


7085 


6A 


2006 Dec 8 


7228 


6B 


2008 Jan 4 


7620 


6A 


18 GHz Observations 






2003 Jul 31 


6002 


6D 


2004 May 6 


6282 


6C 



8640 


12 




8640 


10 




8640, 8900 


13 


1992.9 


8640, 8900 


9 


1992.9 


8640, 8900 


6 


1992.9 


8640, 8900 


8 


1993.6 


8640, 8900 


10 


1993.6 


8640, 9024 


17 


1993.6 


8640, 9024 


9 


1994.4 


8640, 9024 


21 


1994.4 


8640, 9024 


10 


1994.4 


8640, 9024 


12 


1995.7 


8896, 9152 


7 


1995.7 


8640, 9024 


9 


1995.7 


8640, 9024 


11 


1996.7 


8640, 9024 


13 


1996.7 


8896, 9152 


8 


1996.7 


8512, 8896 


7 


1998.0 


8896, 9152 


10 


1998.0 


8512, 9024 


7 


1998.0 


8896, 9152 


12 


1998.9 


8502, 9024 


11 


1998.9 


8512, 8896 


10 


1999.7 


8768, 9152 


11 


1999.7 


8512, 8896 


14 


1999.7 


8512, 8896 


10 


2000.9 


8512, 8896 


11 


2000.9 


8768, 9152 


8 


2001.9 


8512, 8896 


8 


2003.0 


8512, 9024 


9 


2003.0 


8768, 9152 


10 


2003.6 


8512, 8896 


9 


2004.0 


8512, 8896 


9 


2004.0 


8512, 8896 


9 


2004.4 


8512, 8896 


9 


2005.2 


8512, 8896 


9 


2005.5 


8512, 8896 


9 


2006.2 


8512, 8896 


9 


2006.5 


8512, 9024 


8 


2006.9 


8512, 9024 


11 


2008.0 


17345, 19649 


8 




17345, 19649 


7 





t Some datasets have been averaged together to generate the corresponding images in Fig. [T] 
and [2] for the listed epoch. 
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Table 2. Best-fit torus parameters for each epoch, with la confidence intervals. 



Day /(mjy) R (") 6 (°) <5 (%) (°) g (%) X 2 /^ dof§ 



9 GHz Observations 



1786 


3 9+ ' 2 
o.»_ .2 


n cn+0.07 


70 +io 

' u -14 


162+ 19 


188+20 


00 + I6 
°3_26 


1.8 


2115 


1852 


o 7 +0.2 


°- 63 -0.08 


46+19 


108l 6 6 6 2 


57+ 54 
°'-19 


76l!t 


3.7 


1650 


2067 


5 5+ ' 1 
°-°-0.2 


n 70+O.OI 


56l? 8 


10811 2 


143l 3 23 


90+1 

ay -19 


4.3 


3610 


2142 


5 0+ ' 2 
O.U_ 2 


66+ 03 
u - DD -0.02 


63l? 8 


01+101 

81-45 


142115 


45 +54 


17 


2710 


2143 


5 6+ ' 2 
°' D -0.1 


64+ ' 02 
u - D4 -o.oi 


6+ 24 
u.o_ 4 


i+37 
1-1 


09+25 

yz -15 


39+f 


16 


1400 


2313 


6 7+ ' 1 
D -'-o.i 


62+ ' 01 
U ' DZ -0.01 


33l a 


i+33 
1-1 


831?, 


40l 7 


3.6 


2910 


2320 


7 +0 ' 2 

' - u -0.2 


n C7+0.02 
°' 67 -0.03 


31 + 15 


fio+33 
M -48 


78lf 


41+ 9 


4.5 


2970 


2426 


6 7+ ' 1 
D -'-o.i 


690+ ' 009 


53l 2 


36135 1 


75l 9 8 


43l 7 


5.2 


4380 


2550 


6.5+°;* 


n 70 + O.O2 
u - l3 -0.06 


57l 3 25 


01 +85 
yl -23 


133l 9 27 


811- 


7.5 


3000 


2681 


8.4t[- 


672+ ' 005 
u - D ' z -0.005 


46l 2 


™±Z 


78l® 


35l 4 


6.0 


6150 


2685 


O 1+0.1 

8-1-0.1 


665+ ' 007 


52l 2 


1+I6 
1-1 


loelfo 


3711 


5.7 


3264 


3073 


11 l +0 - 2 


o-67il2;22 8 8 


32l| 


52+ 11 


80l| 


4111 


5.6 


2670 


3109 


Q 7+0-4 
y -'-0.3 


0.63312:2? 


16ll« 




751- 


43l 5 3 


17 


1606 


3178 







685+ '" 05 
• DOO -0.005 


44l 2 


n 4+ 3 

u -^-0.3 


94l 4 4 


38.011- 4 


4.0 


3450 


3436 


15.21°;! 





7„- + 0.004 


46.3+?2 


26+ 8 
ZD -14 


95li 


41.911:1 


4.9 


4345 


3485 


15 4+ ' 1 





7n7 +0.004 
• ,u ' -0.004 


SO.!!!" 


4+ 7 

u - -0.29 


102+ 3 


43.3ll;| 


5.2 


4710 


3512 


15 4 +0 1 





7n o + 0.005 
■ ,uo -0.00<i 


5i.ill;2 


1 +25 

U -0.1 


941^ 


4211 


3.3 


2640 


3914 


17 6+ 02 
1, - D -0.2 





504+0.007 


4ll 2 


04+ ' 5 
u - u **_0.04 


112l 4 


3911 


2.4 


1280 


4013 


i9.il° ;l 





?t - 1+0.004 
• ' d -0.004 


45.511:2 


n Q+03 
U - 3 -0.3 


10511 


45.111;= 


3.0 


2838 


4016 


18-712:1 





7.. + O.OO5 
■ '* d -0.004 


49.4+11 


o.il2:i 


IO3I3 


45.7ii : ; 


3.1 


2520 


4220 


20.2+21 





79Q+0. OO-'i 
' ,za -0.0()4 


42.6+ / 9 


6l« 


iool 2 


37 6+ ' 8 


2.2 


2963 


4268 


21.8+°* 





7o K +().()04 
• ,oo -0.004 


40.0li;2 


30l? 5 


107 + 2 


38.012; 9 


7.3 


3870 


4372 


22.912' 2 





797+0.004 
• ,Z ' -0.004 

' ' ' -0.005 


36.6ll; 2 


9 4 +6 
Z4 -10 


11411 


37.5+2'? 


3.6 


3540 


4577 


23.912'j 





39.9l°* 


22+ 9 
zz -19 


10311 


39.4121 


7.0 


3450 


4584 


oc O+0.09 
zo - -0.08 





7 i 7 +0.002 
' -0.002 


41 4+0-> 


4+° 1 
u '^-0.3 


10911 


38.5l2;2 


2.9 


4230 


4966 


90 4+0.05 
■"-0.06 





764 +0.002 
-0.002 


40.ol2; 7 


u " 3 -0.2 


IO8.5+I3 


39.6+2-g 


1.3 


50697 


5011 


oo q+0.07 





7 7K +().()02 
• ' '"-0.002 


43.4+2' 7 


l +0 ' 7 
u - -0.02 


104.9li: 4 


40 1+ ' 6 


1.4 


50539 


5387 


34 1+ 01 
^^-O.l 





7 oq+0.003 
' ' oa -0.()():S 


40 7 +1 '° 
* u -'-0.9 


01+3.4 

u - ul -o.oi 


108.21H 


41 4+0-5 


1.6 


39612 


5748 


41 7+0.O6 

f * J --' -0.08 





011+O.OOI 
01 -0.002 


43.012'j 


10+ 3 ' 4 
u - ±u -0.05 


103.01a! 


40 4+ - 3 
4U - 4 -0.5 


1.3 


39000 


5810 


42 5+ ' 07 





„, .+0.001 
olo -0.001 


42.212^ 


001 + 04 
u.uu±_ 001 


118.112' 9 


43 2+ - 4 
1^- z _0.5 


1.6 


46020 


6003 


46 5 +0 1 
^ D '°-0.08 





01 o+0. 003 
• 01 -0.001 


38.8lJ : ° 


0fi+ 3 ' 7 

u - UD -o.oi 


100-3l2. 04 


00 9+0.5 
oo.z q g 


1.6 


45000 


6129 


52 S+ IU)S 
oz '°-0.07 





000+0.002 


41 0+0.6 


01+0-3 
u - ul _0.01 


107.7l?-2 


42 2 +t) -' 4 


1.3 


46020 


6170 


54 +0 1 





001+0.002 
■ 00 -0.001 


41.4121 


17+ 2 ' 4 

u - ±( -0.14 


108.1+v: 3 


QQ Q + 0.5 
OO.O Q g 


1.5 


43680 


6283 


53.612:2^ 





Q9Q+0-001 
■^"-o.ooi 


39.012J 


16+ 10 

U-4O_0.ll 


i07.3l ;2 


38-7+2-j 


1.5 


44850 


6605 


6i.3i2;l 





011+0.002 

• oyi -0.002 


q 7 7 +0.5 

J '-'-0.6 


19+ 2 - 1 
U - ±M -0.14 


108.2l ;2 


39.0+23 


2.9 


42900 


6693 


62 7 +0 ' 07 
oz - ' -0.08 





oco + 0.001 

• O0 °-o.ooi 


42 2+ ' 4 


11+ ' 8 
u - 11 -0.06 


100.6121 


35-8+21 


1.5 


39000 


6973 


70 7 +0.08 
-0.07 





„oi +0.001 
' 88i -0.001 


43.812'j 


01+ 3 ' 6 

u - ui -o.oi 


116.4+2' 7 


41 8+ ' 4 


1.5 


40365 


7085 


77 1+0.07 
' ' • 1 -0.09 





079+0.001 

'°' -0.002 


38-8l2'.5 


03+ ' 9 
u - uo -o.ooi 


110.512;* 

107.312;* 


39.5+21 


1.5 


29120 


7228 


82 4 +0 ' 10 
° z -^_0.09 





071+0.002 
•°' -0.001 


39.4121 


o-Q7l2:2 4 


39 1+ ' 3 

jm -i_o.2 


1.4 


35685 


7620 


qo c+0.10 
yJ -°_0.10 


o.893i2:221 


41.812-1 


o-o3l2: 00 


103.712J 


38.612; 3 


1.3 


42900 



18 GHz Observations 

6002 25.ll2;2 7 0.80l2;22! 32.ll2;2 0.031213 99.1111 36.5l2;2 0.7 42939 
6282 24.31212 °- 81 lo:oo2 2 6-9_a! 10I9 108.2lj;!j 36.2l°;| 1.0 31200 



^The x 2 values include the statistical errors only, while the uncertainties in the parameters include the 
systematic errors (see text). 

§ Before 2000, all 13 frequency channels in the data were averaged to boost the signal and to reduce the 
computing time. 
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Table 3. Best-fit parameters for the shell fit and shell+2 points fit. 



Fit to shell alone Fit to shell + 2 points 



Day / sh cll R X 2 /" /shell R fl 5RAi iDcc! f 2 SRA 2 SDcc 2 / ratio t di/H§ d 2 /fl§ x 2 /" 
(mJy) (") (mjy) (") (mjy) (") (") (mjy) (") (") 



1786 


3.8 





67 


1.8 


2.5 





97 


1.2 


-0.04 


1852 


3.6 





64 


3.7 


2.6 





57 


0.5 


-0.02 


2067 


5.3 





67 


4.4 


4.0 





81 


1.1 


0.17 


2142 


4.9 





67 


17 


4.1 





76 


0.7 


0.17 


2143 


5.6 





63 


16 


4.2 





67 


1.1 


0.41 


2313 


6.7 





62 


3.6 


5.1 





64 


1.2 


0.36 


2320 


7.0 





67 


4.5 


5.3 





72 


1.2 


0.36 


2426 


6.6 





67 


5.3 


5.3 


o 


76 


1.1 


0.31 


2550 


6.3 





71 


7.5 


4.9 





76 


0.6 


0.77 


2681 


8.4 





66 


6.1 


6.7 





72 


1.3 


0.30 


2685 


8.0 





66 


5.7 


6.5 





70 


1.2 


0.30 


3073 


11 





69 


5.9 


8.8 





70 


1.8 


0.45 


3109 


10 





69 


17 


6.7 





66 


2.2 


0.43 


3178 


12 





69 


4.2 


9.2 





73 


2.0 


0.39 


3436 


15 





69 


5.3 


12 





76 


2.5 


0.35 


3485 


15 





69 


5.4 


13 





76 


2.5 


0.36 


3512 


15 





69 


3.5 


12 





77 


2.6 


0.37 


3914 


17 





72 


2.8 


14 





76 


2.9 


0.43 


4013 


19 





74 


3.6 


15 





80 


3.2 


0.41 


4016 


19 





74 


3.6 


14 





77 


3.5 


0.41 


4220 


20 





73 


2.7 


16 





77 


3.2 


0.41 


4268 


22 





73 


7.9 


17 





80 


3.6 


0.38 


4372 


23 





75 


4.3 


18 





78 


3.6 


0.40 


4577 


24 





76 


7.8 


19 





81 


4.0 


0.42 


4584 


25 





74 


3.6 


20 





79 


4.0 


0.43 


4966 


29 





76 


1.4 


23 





81 


4.8 


0.43 


5011 


33 





77 


1.5 


27 





82 


5.2 


0.43 


5387 


34 





80 


1.7 


27 





83 


5.5 


0.45 


5748 


41 





81 


1.5 


34 





85 


6.5 


0.46 


5810 


42 





82 


1.7 


34 





85 


6.4 


0.47 


6003 


46 





83 


1.8 


38 





85 


7.1 


0.50 


6129 


53 





84 


1.5 


42 





88 


8.6 


0.47 


6170 


54 





85 


1.7 


44 





87 


8.4 


0.50 


6283 


53 





84 


1.8 


44 





87 


7.8 


0.50 


6605 


61 





86 


3.3 


49 





88 


9.5 


0.51 


6693 


62 





87 


1.8 


52 





88 


8.6 


0.54 


6973 


73 





90 


1.9 


60 





92 


10 


0.48 


7085 


77 





89 


2.3 


64 





92 


11 


0.52 
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Note. — /shell and Ft arc the flux density and radius of the thin spherical shell. The point source in the east has a flux 
density j\ and is offset by <5RAi and <5Dcci from the shell center, in R.A and Dec respectively. Similarly, f'2 , 6RA2 and 5Dcc2 
are for the point source in the west. 

'''The flux ratio between the eastern and western lobes is defined as / ra tio — { shc11 j 1 . 

J shell / z * J 2 

§The offset from the shell center to the point source is defined as dj — \/(<5RAi) 2 + (<5Dcc^) 2 . 



